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Abstract: The ANTARES telescope is being built in the Mediterranean Sea. The detector consists of a 
3D array of photomultipliers (PMTs) that detects the Cherenkov light induced by the muons produced 
in neutrino interactions. Other signatures can also be detected. Since the neutrino fluxes from point-like 

■ sources are expected to be small, it is of the utmost importance to take advantage of the ANTARES point- 
ing accuracy (angular resolution better than 0.3 degrees for muon events above 10 TeV) to disentangle 
a possible signal from the unavoidable atmospheric neutrino background. In order to distinguish an ex- 
cess of neutrino events from the background, several searching algorithms have been developed within 
the ANTARES collaboration. In this contribution, the discovery potential and sensitivity to point-like 

' sources of the ANTARES neutrino telescope are presented. 



Introduction 

The ANTARES collaboration [1] has started the 
construction of an underwater neutrino telescope 
in the Mediterranean Sea at a depth of 2475 m. 
Seven of the twelve lines of the detector have been 
deployed so far (May 2007) and the whole de- 
tector will be commissioned in early 2008. Each 
line is equipped with 75 10" photomultiplier tubes 
(PMTs) joined in triplets making 25 floors along 
the line. The lines are kept vertical by means of 
a buoy located at their upper end. The mean dis- 
tance among lines is about 65 m and the instru- 
mented length starts at 100 m over the seabed and 
covers about 350 m. The angular resolution for 
the ANTARES telescope at high energies (above 
10 TeV) is better than 0.3°. At lower energies, the 
angular resolution is dominated by the angle be- 
tween the muon track and the original neutrino di- 
rection. The angular resolution together with the 
effective area determine the performance of a neu- 
trino telescope. The good pointing accuracy and 
large effective area at high energies (0.06 km^ at 
100 TeV) enables us to search for point-like neu- 
trino sources or, in case no hint of neutrinos source 



is found, to set restrictive upper limits to neutrino 
fluxes. 

Search for point-like sources 

In the search for point-like steady sources, the 
whole visible sky from the ANTARES location is 
surveyed making no assumption about time cor- 
relation with observations of any other physical 
phenomenon, like GRBs, flares and other tran- 
sient sources. In this contribution, we will intro- 
duce a new method for the search for point-like 
sources. This method is based on the Expectation- 
Maximisation (EM) algorithm which is widely 
used as a likelihood maximisation algorithm for 
clustering analysis. 

The EM algorithm 

The Expectation-Maximisation [2] algorithm is a 
general approach to maximum-likelihood estima- 
tions for finite mixture model problems. For this 
method a parametrisation of both the signal and 
neutrino background density distributions is re- 
quired. The main assumption is that sources sig- 
nals are supposed to follow Gaussian distributions 
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which is reasonable assumption in our case. The 
background distribution is inferred from the Monte 
Carlo data sample, but it could be obtained directly 
from the real data by scrambhng the right ascen- 
sion coordinate of the measured events. No energy 
information or further performances of the detec- 
tor are used in this case. The parameters are de- 
termined maximising the likelihood using the EM 
algorithm. The EM method assumes the existence 
of missing or hidden information, in our case the 
knowledge whether a neutrino event comes from a 
given source or it is produced by the atmospheric 
neutrino background. Hence, the real observed 
data can be understood as an incomplete sample, 
so that, adding a new vector z we can build the 
complete data sample, where z is a class indica- 
tor vector that tells whether an event comes from 
a source or not. Maximisation of likelihoods an- 
alytically intractable can be easily accomplished 
by means of this methodology. The EM general 
method follows an iterative procedure where each 
iteration has two steps: 

1. E-step: The expected value of the complete 
data log-Ukelihood, conditional to the ob- 
served data, is computed for a given set of 
parameters 

2. M-step: Find * = ijr(™+i) that maximises 
the expected value. This maximisation will 
lead to the maximisation of the desirable 
log-likeUhood of the complete data sample 

The model selection criteria is performed by means 
of the Bayesian Information Criterion or BIC. 
This BIC value is an approximation of the inte- 
grated likelihood when the number of events is 
high enough. The BIC value can be used as a test 
statistic and the discovery potential or discovery 
power of a point-Uke source to be detected can be 
inferred using the hypothesis testing theory. This 
discovery power wiU be described in section 3. 

Other searching algorithms 

In ANTARES, other clustering algorithms have 
been developed in order to exploit the potential of 
the detector in the search for point-like sources. 
The first approach was a binning method used to 
look for a cluster of events in different bins of a 



grid in which the sky is divided [3]. The signifi- 
cance of clusters is estimated comparing with the 
distribution of atmospheric neutrino events which 
is uniform in right ascension. These methods rely 
on the knowledge of the angular resolution of the 
detector in order to build the optimum grid. Nev- 
ertheless, no hypothesis is made on the neutrino 
source apart from the neutrino emission according 
to a power law with spectral index of 2. 
The first algorithm that did not rely on a binning 
approach was a Likelihood Ratio method (LR) [4]. 
This method operates by testing the data compat- 
ibility with two hypotheses: one, often called Hq, 
is the null hypothesis and the other one called Hi. 
The null hypothesis is assumed to be the presence 
of background only, and Hi considers the existence 
of a neutrino point source in addition to the atmo- 
spheric background. The parametrisation of the 
density functions are based on the obtained infor- 
mation from the simulated expected performances 
of the detector. 

Foreseen results for the search for 
point-Uke sources in ANTARES 

In this section, we will review the results of the 
point-like source search in ANTARES based on 
a blind search over one year of data taking. A 
blind search is called when no assumption about 
the source is done beyond the educated assump- 
tion of a power law neutrino emission from the 
sources. This type of search differs from the fixed 
point searches, where a set of candidate sources are 
studied. The most direct way of presenting the re- 
sults in the searching analysis is by means of the 
discovery power, also called discovery potential. 
This parameter accounts for the percentage of suc- 
cess of discovering a point-like source over the at- 
mospheric neutrino background. Therefore given a 
number of events per year observed at the detector 
from a source, the power of the method is the prob- 
ability of detecting such source with a given con- 
fidence level. This definition enables the construc- 
tion of plots as the one shown in figure 1 . This fig- 
ure shows the discovery potential for a confidence 
level of 5(7 as a function of the mean number of ob- 
served events from a source located at the vicinities 
of (5 = —80°. The three main searching algorithms 
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- EM algorithm (S = -80°) 
LR algorithm (S = -81°) 

- Binning method (5= -73.4°) 
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Figure 1: Discovery power (probability of detec- 
tion) for a source located at low declinations as a 
function of the averaged detected events for three 
methods developed in ANTARES for 5a confi- 
dence level. 

are shown. It is important to mention that the re- 
sults yielded by the different analysis correspond to 
different Monte Carlo simulations, muon track re- 
construction strategies, and therefore different per- 
formances of the detector. Hence a direct compar- 
ison of the results obtained is not completely fair. 
Nonetheless, we can observe that unbinned tech- 
niques result in a sound improvement with respect 
to the classical binning approach. This result mo- 
tivates the use of unbinned techniques in the point- 
like analysis of real data. Moreover, in this work, 
the energy estimate information is not used in the 
EM algorithm which makes it specially appealing 
for the initial years of operation in a neutrino tele- 
scope where the information from the muon energy 
estimate does not constitute a solid piece of knowl- 
edge. 

The discovery potential depends on the source lo- 
cation in the sky, specially it depends only on the 
source declination since the atmospheric neutrino 
background is right ascension independent. In or- 
der to show more clearly this dependence, it is 
commonly presented as the required number of ob- 
served events from a source to claim its existence 
in 50% of equivalent experiments as a function of 
the declination. Figure 2 shows this magnitude for 
two different confidence levels using the EM al- 



Figure 2: Mean number of observed events re- 
quired to yield a discovery power of 50% for dif- 
ferent confidence levels as a function of declination 
using the EM algorithm. 

gorithm. The required number of events is larger 
at lower declinations since the number of atmo- 
spheric neutrinos per solid angle is larger at low 
declinations. In order to present the results in terms 
of a neutrino/muon flux, we can divide the required 
number of events by the exposure of the detector in 
one year to that declination. Figure 3 shows the 
required integrated neutrino flux needed to claim 
the existence of a neutrino source with a probabil- 
ity of 50% as a function of the declination. The re- 
quired flux increases with declination mainly due 
to the visibility factor. From the ANTARES loca- 
tion, lowest declinations are visible 100% of the 
time, whereas high declinations have a small vis- 
ibility, so that higher fluxes are needed to make a 
discovery. 

Sensitivity 

If no excess of events is observed in the sky and 
therefore, no evidence of a neutrino point-like 
source is claimed, we can set an upper limit on 
the expected flux from a neutrino emitter This 
upper limit can be provided in terms of a full sky 
blind search. Nonetheless, most commonly it is 
presented in terms of a point to point upper limit. 
Moreover, the concept of upper limit is only ap- 
plied when the experiment is currently running, an 
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Figure 3: Energy integrated (above 10 GeV) neu- 
trino flux needed to claim the existence of a source 
with probability of 50% for two confidence levels 
and the EM algorithm for a E^'^ spectrum. 

averaged upper limit or sensitivity where the aver- 
age upper limit is obtained by an ensemble of 10^ 
experiments. Figure 4 shows the expected sensi- 
tivity for ANTARES for one year of data taking 
compared to the published results from other ex- 
periments. MACRO [5] results are computed after 
alive-time of 6.3 years. AMANDA upper limits 
are taken from [6, 7]. The projected sensitivity of 
IceCube averaged over all Northern declinations is 
shown [8]. The required flux for a 5ct discovery at 
50% for ANTARES in one year is also indicated. 

Summary and conclusions 

In this contribution we have reviewed the expected 
performance of the ANTARES neutrino tele- 
scope regarding the search for point-like neutrino 
sources. Several searching algorithms have been 
devised in ANTARES. Among them, unbinned 
techniques turn out to be more efficient than the 
standardised binning approach. In addition, the 
method based on the EM algorithm presents very 
good results without the trade-off of higher de- 
pendence on the estimated detector performances. 
Discovery potentials in terms of number of events 
and required neutrino flux to claim the existence of 
a source at 50% of probability has been presented. 
The comparison with other experiments in term of 



Figure 4: Sensitivity of ANTARES in one year 
compared to other experiments. See text for fur- 
ther details. 

the sensitivity has also been presented. 
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